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Abstract. We discuss XMM-Newton and BeppoSAX observations of MKN 231, the lowest-redshift Broad 
Absorption Line (BAL) QSO known so far and one of the best-studied Ultra Luminous Infrared Galaxies. By 
combining the XMM-A^euiton spectral resolution and the high-energy bandpass of BeppoSAX we have been able 
to study in more detail than previously possible its 0.2-50 keV spectral properties. The 5eppoSAX PDS data 
unveiled, for the first time, a highly absorbed (TVh 2 x 10^'*cm~^) power-law component. We find that: a) below 
10 keV we are seeing only reprocessed radiation through reflection and/or scattering; b) the intrinsic 2-10 keV 
luminosity of MKN 231 is l^o's x 10** ergs s~^, i.e. more than an order of magnitude greater than previous 
measurements; c) the starburst activity significantly contributes to the soft {E < 2 keV) X-ray emission. 

Key words. Galaxies: active. Galaxies: starburst. X-rays: galaxies, Galaxies: individual: Markarian 231 



^ 1. INTRODUCTION starburst and, in particular, radio observations of the in- 

' ner kiloparsec regions suggest that a nuclear ring of active 

X ■ In the local (z < 0.05) Universe Markarian 231 (z ^ star formation is proceeding at a rate of - 220 Mq yr'^ 

0.042, Lboi = 6 X lO^^Lg for Hq = 50 km s^^ Mpc^^) is (Taylor et al. 1999)^ Indication of a strong starburst ac- 

■ - - ■ one of the most luminous and best studied Ultra Luminous tivity, which can account for up to 70% of the Lpm of 

Infrared Galaxies (i.e. sources with Lm > IO^^Lq; here- mKN 231 (corresponding to an SFR ~ 470±50 Mq yr"!), 

after ULIRGs). This object is optically classified as a comes also from the modelling of the radio to optical spec- 

Broad Absorption Lme (BAL) QSO (Smith et al. 1995). tral energy distribution (Farrah et al. 2003). This latter 

This class of sources, which represent - 10% of the opti- estimate is in good agreement with the AGN fraction de- 

cally selected QSOs, are known to be weak X-ray sources, ^ived by Lonsdale et al. (2003) using VLBI observations 

their faintne ss bei ng attributed t o intr insic absorption at 18 cm. As discussed by several authors the combination 

(Green et al.|2001J Gallagher et al.|2002I). of starburst and luminous AGN makes MKN 231 one of 

Although the nature of the primary energy source the best examples of the transition from starburst to AGN 

(AGN vs. starburst activity) for the high luminosity of according to the scheme outlined in Sanders et al. H199f)|l . 

ULIRGs is still a debated issue (see Ptak et al. 2003, The strong circumnuclear starburst activity together 

Franceschini et al. 2003 and references therein) , in the case with the presence of high intrinsic absorption, probably as- 

of MKN 231 optical and near infrared observations seem sociated with the BAL outflows (Green et al. 2001), makes 
to suggest that a significant contribution to the infrared 

luminosity (LpiR c:^ 3 x IO ^^Lq) could be ascrib ed to AGN i in this respect it is worth noting that the starforming rate 

activity (Goldader et al. 119951 Krabbe et al. 1199711 . On in M82, the prototype local starburst galaxy, is about one order 

the other hand this galaxy is also undergoing an energetic of magnitude lower (Grimm et al. 12003^ 
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the X-ray continuum of this source particularly complex. 
Many different spectral components should be taken into 
account in the spectral modeling i.e. the thermal and the 
binary emission associated with the products of the star- 
burst as well as the scattered/reflected and intrinsic com- 
ponent(s) associated with the AGN. 

Previous X-ray observations by ROSAT, ASCA 
(Turner 1999, Iwasawa 1999; Maloney & Reynolds 2000, 
hereafter MROO) and more recently Chandra (Gallagher 
et al. 2002, hereafter G02; Ptak et al. 2003), revealed 
the presence of: a) extended soft X-ray emission of ther- 
mal origin, likely associated with the circumnuclear star- 
burst; b) a hard and flat (photon index F ~ 0.7) non- 
thermal component having an observed X-ray luminosity 
L(2-io) ^ lO^'^ergs s~^ and c) a Fe Ka emission line hav- 
ing an observed equivalent width of ~ 300 eV (the line 
was measured by ASCA but is not detected by Chandra) . 

To explain the low observed luminosity (which is 
at least two orders of magnitude lower than that ex- 
pected for an AGN of similar bolometric luminosity), 
the flat photon index and the EW of the Fe Ka line, 
MROO suggested that about 3/4 of the observed 2-10 
keV spectrum of MKN 231 consists of a scattered 
power-law component plus a reflected component (no 
direct AGN component is seen below ~10 keV). Both 
reflected and scattered components are absorbed by 
an intrinsic column density of A^h 3 x 10^^ cm~^. 
On the contrary G02, using Chandra data (which in 
addition make it possible to disentangle the nuclear from 
the extended component), support a model in which a 
small, Compton-thick absorber blocks the direct X-rays. 
Only indirect scattered emission from multiple lines 
of sight can reach the observer in the Chandra energy 
range; in particular they proposed a triple scattering 
model composed of three separate absorbed power laws 
having a photon index fixed to 2.1 and absorbing column 
densities of ~ 1.2 x lO^^cni-^ A^h ^ 3.2 x lO^^cm'^ 
and A'h — 5.9 x lO^^cm^^, respectively. According to 
the modeling of MROO and G02, the X-ray emission 
below 10 keV is only reprocessed radiation, thus they 
could evaluate the intrinsic luminosity of MKN 231 only 
through indirect arguments. 

Here the results from XMM-Newton and BeppoSAX 
observations of this object are presented. The main result 
from these data is that a highly absorbed (A^h ~ 2 x 10^^ 
cm~^) AGN component having an intrinsic 2-10 keV lumi- 
nosity of 1^0.5 X lO'*'* ergs s~^ has been detected. This is 
the first direct measurement of the intrinsic power 
of the AGN hosted by MKN 231. This paper is or- 
ganized as follows: observation and data analysis are de- 
scribed in Sect. 2; timing and spectral analysis are re- 
ported in Sect. 3 and 4. Results and discussion are pre- 
sented in Sect. 5. In this paper we assume Hq = 50 km 
s-i Mpc"^ and qa = 0.5. 



2. Observations and Data reduction 

MKN 231 is part of an XMM-Newton survey of 10 nearby 
(z < 0.2) ULIRGs for which high-quality mid-IR and 
optical spectroscopic data are available. Results on this 
project are presented and discussed in Braito et al. (2002) 
and in Franceschini et al. (2003). MKN 231 also belongs 
to a small sample of ULIRGs for which we have obtained 
BeppoSAX measurements. 

2.1. XMM-Newton data 

MKN 231 was targeted for 20 ks by XMM-Newton on July 
6th, 2001. The EPIC cameras (European Photon Imaging 
Camera: IStriider et al. 20011 and ITurner et al.-20?rni were 
operating in full-frame mode with the thin filter applied. 
Event files produced from the standard pipeline processing 
have been analyzed and filtered (using version 5.4 of the 
Science Analysis Software, SAS, and the latest calibration 
files released by the EPIC team) from high background 
time intervals and only events corresponding to pattern 0- 
12 for MOS and pattern 0-4 for pn have been used (see the 
XMM-Newton Users' Handbook, Ehlc et al. 2001|l . The 
net exposure times, after data cleaning, are ~ 16.0 ksec, 
- 19.9 ksec and - 19.8 ksec for pn, MOSl and M0S2 
respectively. 

In the low-energy domain (E < 2 keV), the X-ray 
emission of MKN 231 appears to be extended^, with a 
roughly circular symmetry. The presence of extended soft 
X-ray emission was first noted using ROSAT observations, 
(Turner T9W) and it has been confirmed by the recent 
Chandra observation (G02). The Chandra data show that 
the galaxy emission (0.35-2 keV range) extends over a ra- 
dius of ~ 25 arcsec. The hard (2-10 keV) X-ray brightness 
profile is comparable to the XMM-Newton PSF, whose 
HEW at ~ 8 keV is ~ 14". 

Since in this paper we want to address the X-ray spec- 
tral properties of the nuclear and circumnuclear X-ray 
emission, source counts have been extracted from a cir- 
cular region of 20 arcsec radius'^, positionally coincident 
with the core of the hard X-ray emission. The X-ray po- 
sition is also coincident (within the XMM-A^ewton posi- 
tional accuracy ^ few arcsec) with the core of the optical 
and radio emission. The XMM-Newton PSF does not al- 
low us to select a smaller region and thus, although the 
extraction radius has been chosen to minimize the effect 
of the extended soft emission, in the derived spectrum we 
expect imprinted signatures of the strong circumnuclear 
star burst activity of MKN 231. 

The presence of other sources around MKN 231 pre- 
vented the selection of a nearby large enough area to accu- 

A detailed spectral and imaging analysis of the extended 
component will be reported elsewhere. 

The X-ray source detected by G02 using Chandra data at 
3" apart from the nucleus is thus included in the counts ex- 
traction region. However, as its flux is two orders of magnitude 
lower than that of the nuclear emission of MKN 231, it does 
not affect the results reported here. 
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mulate background counts with good statistics; therefore 
background counts have been accumulated using the back- 
ground event files distributed by the EPIC team. Response 
matrices (that include the correction for the effective area) 
have been generated using the SAS tasks arfgen and rm- 
fgen. The net count rate in the 0.2-10 keV energy range 
is (11 ± 1) X 10~2 cts/s, (3.77 ± 0.14) x 10"^ cts/s and 
(3.54 ± 0.14) X 10-2 cts/s for pn, MOSl and M0S2, re- 
spectively. 

2.2. BeppoSAX data 

MKN 231 was observed by BeppoSAX (Boella et al., 
nWf^ from December 29th, 2001 to January 1st, 2002. 
In this paper we use data collected from the Medium 
Energy Concentrator Spectrometer (MECS) and from the 
Phoswich Detector System (PDS). We do not use data 
from the Low Energy Concentrator Spectrometer since 
they are of significantly lower quality when compared with 
the low energy XMM-Newton data. 

Standard data reduction was performed 
using the s oftware package SAXDAS (see 
jhttp : // www . sdc . asi . it / software and the Cookbook; 
Fiore, Guainazzi & Grandi 1999) for the MECS datasets 
and XAS (Chiappetti & Dal Fiume 1997) for the PDS 
data. Spectral fits were performed using public response 
matrices as from the 1998 November issue. BeppoSAX 
pointed at MKN 231 for an effective exposure time of 144 
ksec and 76 ksec for MECS and PDS, respectively. To 
maximize statistics and signal-to-noise ratio, the MECS 
source counts were extracted from a circular region of 2 
arcmin radius. The background spectrum was taken from 
blank sky fields at the position of the source with the 
same extraction radius. As suggested in the BeppoSAX 
Cookbook we have considered only the MECS data in 
the 1.8-10 keV energy range. Spectral channels corre- 
sponding to energies in the range 15-60 keV were used 
for the PDS. MECS source counts have been rebinned 
proportionally to the instrumental resolution and so as 
to have always more than 20 counts per channel; the 
PDS data instead have been rebinned to maximize the 
signal to noise ratio in each energy bin. The net count 
rates are (5.68 ± 0.23) x lO'^ cts/s in the MECS and 
(6.58 ± 2.16) X 10-2 cts/s in the PDS. 

3. X-ray variability 

The Chandra data of MKN 231 discussed by G02 indicate 
significant variability with a count rate decrease of ^ 45% 
on a 20 ks time-scale at energies above 2 keV. This vari- 
ability coupled with the 188 eV upper limit on the EW of 
the Fe Ka emission line was used by G02 to argue against 
the reflection model proposed by MROO. 

Thus we first investigated the variability properties of 
MKN 231 using the BeppoSAX data. Although BeppoSAX 
has a significantly worse spatial resolution than Chandra, 
the long exposure time (144 ks) of this pointing is ideal 
to look for variability on similar time scales as those in- 



vestigated by G02. Using the BeppoSAX MECS counts in 
the 2-10 keV energy range we do not find evidence of 
variability: the test, applied to the binned light curve 
(AT — 5700 sec), gives a probability of 50% that the data 
are consistent with a constant flux during the BeppoSAX 
observation. It is worth noting that this result does not ex- 
clude variability on different time scales and/or irregular 
variability. 

Finally, the XMM-Newton and BeppoSAX observa- 
tions are not simultaneous but the MECS and XMM- 
Newton 2-10 keV fluxes, derived assuming a simple power- 
law (PL) model, are comparable within the relative uncer- 
tainties; thus a statistically significant flux difference is not 
found between the two data sets. 

4. Spectral analysis 

The XMM-Newton and BeppoSAX data have been 
analyzed using standard software packages (FTOOLS 
5.2, XSPEC 11.2). All models discussed here have been 
filtered for the Galactic absorbing column density along 
the hne of sight (A^h = 1-3 x lO^o cm'^ Dickey & 
Lockman, 1990). Unless otherwise stated, errors will be 
given at the 90% confidence level for one interesting 
parameter (Ax^ = 2.71) and the line(s) energy will be 
given in the source rest frame. The thermal emission 
has been described with the MEKAL model in XSPEC 
(Mewe et al. 1985). Abundances were kept solar; keeping 
them free did not substantially change the results. Since 
MKN 231 does not show significant flux variation we have 
first derived the cross-calibration between MOS, pn and 
MECS and then, by assuming that the cross-calibration 
factor between MECS and PDS ranges from 0.7 to 0.95 
(see Fiore, Guainazzi & Grandi 1999), we have tied 
together the XMM-Newton data and the PDS data. 

From previous X-ray observations, performed with 
ASCA and Chandra, we know that the broad-band (0.2-10 
keV) spectrum of MKN 231 appears to be more complex 
than a single power-law model. In particular the Chandra 
data presented by G02 show that the circumnuclear star- 
burst thermal emission, which extends up to ~ 25" from 
the nucleus, dominates the spectrum at i? < 2 keV. Due 
to the PSF of XMM-Newton and BeppoSAX detectors the 
extracted counts in the low energy domain are therefore 
"contaminated" by the circumnuclear starburst. In con- 
trast the hard (E > 2 keV) X-ray emission should be less 
affected by the starburst activity and more suitable for in- 
vestigating AGN related emission(s). For these reasons we 
have decided to focus our attention first on the 2-50 keV 
spectrum, and for simplicity (since we have enough statis- 
tics) to use only the pn and the PDS data. In a second 
step we have also considered the 0.5-2 keV source counts 

^ The BeppoSAX PSF does not allow us to extract a nuclear 
light curve, as done by G02 with the Chandra data, but we 
already know that the soft starburst contribution is not signif- 
icant above 2 keV. 
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as well as the data from all the other detectors (MECS, 
MOSl and M0S2). 

4.1. Fitting the Hard (2-50 keV) X-ray Energy Band: 
pn + PDS data 

In agreement with previous analysis (Iwasasawa 119991 
Turner 1999, MROO and G02) when we fit the 2-10 keV 
spectrum of MKN 231 with a single unabsorbed or ab- 
sorbed PL model, we find a photon index (F = 0.67 ±0.11 
and r = 0.92t°J^, respectively) which is unusually low 
if compared with the typical photon index of radio-quiet 
QSO (r ranges from 1.7 to 2.3, Reeves & Turner IMIHl 
George et al. 2000). Furthermore the fit leaves positive 
residuals around 6.4 keV, strongly suggesting the presence 
of the iron emission line(s). The best fit energy position of 
this line is 6.40tg keV with an observed EW = 300±130 
eV, in agreement with the ASCA data discussed by MROO 
(in this first modelling the Fe line was constrained to 
be unresolved). The flatness of the observed X-ray spec- 
trum is suggested also by the high energy (E = 10-50 
kcV) PDS data: in fact if we constrain the photon index 
to lie in the range 1.7-2.3, the PDS data are underesti- 
mated by at least a factor of 4. These results suggests: 
a) the presence of a heavily absorbed component and b) 
that below 10 keV we are probably seeing only the re- 
processed radiation (through reflection and/or scattering, 
see Matt et al. 1991). In the following we investigate in 
turn two extreme possibilities: a reflection-dominated and 
a scattering-dominated scenario. 

4.1.1. The reflection-dominated scenario 

In the reflection-dominated scenario the bulk of the 2-10 
keV emission derives from Compton reflection from op- 
tically thick material. When we consider reflection from 
neutral matter a strong (i.e. EW as high as a few keV; 
Matt et al. 1996) Fe emission line is also expected, this 
line being produced by the same medium which reflects 
the primary continuum. The lack of a strong Fe emission 
line seems then to argue against a reflection-dominated 
scenario. However, the intensity and the energy centroid of 
the Fe line depends on the ionization state of the medium 
(Matt et al. 1993, Ross et al. 1999); in particular, the in- 
tensity decreases, compared to reflection from a neutral 
medium, when the ionization parameter (^) is between 
100 and 300, due to the resonant trapping opacity (Zycki 
& Czernv fTOMil . 

a) Neutral medium 

A pure reflection-dominated continuum, where the 
reflecting medium is neutral, cannot reproduce the 2-50 
keV spectral properties of MKN 231, leaving strong 
positive residuals above 7 keV. Furthermore this model 
would require a Fe emission line at 6.4 keV stronger 
than that measured. As previously noted by MROO, a 



prominent scattered component is needed to reproduce 
the observed continuum and to dilute the Fe line. We then 
tested the pure reflection plus scattering model using a 
PL plus a pure reflected continuum (PEXRAV^ model in 
XSPEC; Magdziarz & Zdziarski 1995) plus an unresolved 
gaussian line. The best fit photon index (F ^ 1.3) is 
still low if compared with the typical photon index of 
radio-quiet QSOs. Constraining the photon index to lie 
in the range 1.7-2.3 we cannot find an acceptable fit. The 
results are summarized in Table 1, model A and model 
B. These models do not consistently explain the low 
observed EW of the line and the flatness of the spectrum 
for an intrinsic standard AGN emission. 



b) Ionized medium 

An alternative explanation for the observed low EW 
of the Fe line invokes the ionization state of the reflecting 
medium. 

Since the reflected continuum from mildly ionized ma- 
terial predicts the presence of structured features in the 
low-energy domain (see e.g. Ross & Fabian. 1993) we are 
forced to consider also the pn counts below 2 keV. On the 
other hand, below 2 keV we also have to take into account 
the contribution from the circumnuclear starburst. 

Thus, we have first tested a 3-component model com- 
posed of: a) thermal emission to take into account the 
starburst contribution at low energies; b) a pure reflected 
component from mildly ionized material (PEXRIV^ model 
in XSPEC; Magdziarz & Zdziarski 1995) and c) a gaussian 
line. 

This model gives a poor representation of the broad 
band (pn + PDS data) properties of MKN 231 {x^/v = 
102.6/65; see Table 1 model C), leaving strong positive 
residuals at energies above 7 keV and underestimating 
by a factor 5 the counts in the PDS energy range; all 
these features strongly suggest the presence of a heavily 
absorbed (A'h ^ 10^'* cm~^) component. We have thus 
added a highly absorbed PL component (described by the 
PLCABS^ model in XSPEC; Yaqoob 1997) whose photon 
index has been tied to the photon index of the reflected 
component from ionized material. 

This model gives a better broad band description of 
the spectrum of MKN 231 (see Table 1 model D), but it 
requires a very steep intrinsic photon index (F = 2.95 ± 
0.07). 

^ Given the present statistics above 20 keV we fixed cos(i) — 
45, Z = Zq and Ec — 200 keV and allowed only the photon 
index and the normalization to vary. 

® Given the present statistics above 20 keV we fixed cos(i) — 
45, Z = Zq and Ec = 200 keV and allowed only the photon 
index and the normalization to vary. The ionization parameter 
^ was constrained to lie in the range 100-300. 

Given that the intrinsic A'^h is expected to be high (~ 
10^*cm~^) we have used this model, which takes into account 
also the Compton scattering. 
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Finally, given that MKN 231 has a strong circumnu- 
clear starburst (see Introduction), and given the size of 
the source counts extraction radius, we expect a contri- 
bution from X-ray binaries, and in particular from High- 
Mass X-ray binaries (Franceschini et al., 2003; Persic et 
al., in preparation). This latter component has been de- 
scribed assuming an absorbed cutoff-power law model of 
the form exp^'^'^^'^'^ with photon index F ~ 1.1 and 
cutoff energy kT - 10 keV (cfr. Persic & Raphaeh 2002). 

This final model gives a good description of the X-ray 
spectra of MKN 231. In Fig. 1 (left panel) we show the 
data and the residuals while the best fit spectral parame- 
ters are reported in Table 2 (Model A). 

The line energy position is 6.38 ± 0.19 at 90% con- 
fidence level, in agreement with the value (6.5 keV) ex- 
pected for this range of ionization of the mirror. The EW 
of the line is ~ 200 eV with respect to the refiected compo- 
nent, which is consistent with being produced by a mildly 
ionized mirror (Ross et al. 1999). By allowing the line 
width to vary we obtain cr ^ 0.1 keV, consistent with the 
XMM-Newton spectral resolution. Thus at the spectral 
resolution of XMM-Newton the detected Fe line is unre- 
solved. No statistical difference for the best fit parameters 
reported in Table 2 is found if the line width is allowed to 
vary. 

In conclusion, a pure reflection-dominated continuum 
from neutral material cannot consistently explain the low 
Fe emission line EW and the flatness of the spectrum. A vi- 
able solution is the pure reflection continuum from slightly 
ionized material combined with a highly absorbed compo- 
nent. Adding to this "AGN baseline model" the spectral 
components which are likely associated to the starburst 
(a soft thermal component and a PL component associ- 
ated to the High Mass X-ray Binaries) we found a value of 
~ 2.5 for the AGN intrinsic photon index. This steep in- 
trinsic photon index is not so striking if we consider that 
MKN 231 shares many spectral characteristic of NLSls 
(i.e. weak [O777] , relatively narrow Balmer emission lines 
and strong optical Fe// emission; Mathur 2000; Brandt 
& Gallagher I20()()|l . According to this model, hereafter 
called the "AGN reflection-dominated" model, the intrin- 
sic 2-10 keV luminosity of the AGN powering MKN 231 
is L(2_io) 2 X lO'^'' ergs s~^, in agreement with the value 
predicted by MROO based on ASCA data. 

4.1.2. The scattering-dominated scenario 

Another possible explanation for the hard MKN 231 X-ray 
continuum flatness and the lack of a strong Fe emission 
line is a scattering-dominated scenario in which we see the 
AGN mainly through scattered emission (see G02 model- 
ing of the Chandra data). In this context we assume that 
the Fe emission line is produced by transmission and is 
diluted from the scattered component. 

To test this scenario we have used a 3-component 
model composed of: a) a highly absorbed PL represent- 
ing the intrinsic AGN emission transmitted through the 



Compton-thick screen; b) a PL having the same pho- 
ton index as the transmitted component, representing 
the scattered AGN emission and c) a gaussian line. This 
model gives a fairly good description of the data {y^ /v — 
20.47/27) and, as shown in Fig. 1 (right panel), no strong 
residuals are present. 

The intrinsic AGN emission is described by a PL 
model having a photon index of F = 1.69iQ'7Q which 
is filtered by an intrinsic absorption column density of 
A^H = l-4to'5 X lO^^cm^^. The energy position of the 
Gaussian line best fit value \s E = 6.39^0 keV, which is 
consistent with the neutral Fe-Ka line. Allowing the line 
width to vary we obtain a = 0.2 ± 0.18 keV (we note how- 
ever that if we allow the line width to vary the fit does not 
improve significantly). The Fe line EW is ~ 260 eV with 
respect to the observed continuum and ^ 1.2 keV with 
respect to the absorbed PL. This latter value is consistent 
with that expected to be produced by transmission in the 
same medium that absorbs the primary nuclear emission 
(Leahy & Creighton 1993). The scattered component is 
also absorbed by an intrinsic absorption column density 
of Nu = X W^cm-^. 

In the 2-10 keV band the observed flux is 7 x 10~^^ 
ergs cm~^ s~^ and the intrinsic 2-10 keV luminosity, de- 
absorbed for the intrinsic A^h, is L(2-io) = 4 x 10^'^ ergs 

In summary this model provides another possible ex- 
planation for the flat MKN 231 observed 2-50 keV con- 
tinuum and its line properties. 

Furthermore, although not well constrained, the best 
fit photon index is now fully consistent with the range usu- 
ally found in AGN (F ~ 1.7-2.3). In the following we will 
refer to this model as the "AGN scattering-dominated" 
model. 

4.2. The broad band (0.5-50 keV) X-ray spectral 
properties of MKN 231 

In the previous section, using the pn and the PDS data, we 
found that two competitive models (refiection-dominated 
and scattering-dominated) can describe the MKN 231 
AGN emission. In this section these two models will be 
used as starting points to study the broad band (0.5-50 
keV) X-ray spectral properties of MKN 231. Low energy 
data {E < 2 keV) and aU the detectors MOSl, M0S2, 
pn, MFCS together with the PDS data will be used. 

In the ionized refiection model case we found that 
the model derived for the pn and the PDS data (which 
already include the spectral components expected from 
the starburst, see Sect. 4.1.1) is able to reproduce all data 
with no strong differences for the previously determined 
best-fit parameters. In Fig. 2 (left panels) we show the 
model and the residuals for all the detectors, while the 
best fit spectral parameters are reported in Table 2 
(Model A). Within the framework of this modelling the 
intrinsic luminosity of the AGN powering MKN 231 is: 
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Fig. 1. Left panel: 0.2-50 keV pn (filled triangles) plus PDS (empty squares) data and ratio between data and the 
best fit reflection model (model A in Table 2). Right panel: 2-50 keV pn (filled triangles) plus PDS (empty squares) 
data and ratio between data and the scattering dominated model discussed in Sect. 4.1.2 



L(o.5-2) = 3.3 X lO'^^ ergs s ^ and L(2_io) = 2 x lO'^^ ergs 
s-i (see Table 3). 

In the scattering-dominated scenario case the residuals 
below 2 keV are indicative of thermal emission. Indeed, 
two thermal components are needed to reproduce the soft 
X-ray emission of MKN 231 {kT 



0.35+0 ;;^ keV and 



kT ~ 0.91I^;ig keV). For the HMXBs emission we have 
added an absorbed cutoff PL component, following what 
already done in Sect. 4.1.1; moreover, this component is 
needed to reproduce the continuum detected in the 1-3 
keV range. The best-fit spectral parameters derived with 
this broad-band, multi-component modelling are reported 
in Table 2 (Model B; xV^^ 212.9/187). Note that the 
best-fit photon index for the AGN component is now 
F = 1.83t";gg. The derived intrinsic 2 10 kcV luminosity 
for the AGN is ^(2-10) — 4.5 x lO^'^ ergs s^^ A summary 
view of the model is reported in Fig. 2 (right panel). 



5. Discussion 

5.1. The AGN component in MKN 231 

The AGN spectral behavior can be described with two 
models which differ mainly at E ^ 8 keV. Both these 
models require a heavily absorbed component to account 
for the X-ray emission above 10 keV detected with the 
PDS. The photon indices derived with these models, 
although not well constrained, are consistent with the 
typical values of AGNs. Furthermore the PDS data have 
proven, for the first time, the presence of a thick column 
density screen (TVh ^ 2 x 10^** cm~^) which blocks the 
primary X-ray emission. The 2-10 keV observed spectrum 
could be the result of reflection (reflection-dominated 
scenario) or of scattering (scattering-dominated scenario) 



which emerges from the Compton-thick screen. Clearly 
the real physical situation could be an intermediate 
one where we see both components. Since the com- 
plexity of the 2-10 keV continuum prevents us from 
distinguishing, on a statistical basis, between these 
two models, we have checked if the physical parame- 
ters derived in both scenarios are in agreement with 
what is known about ULIRGs, starburst galaxies and 
in particular with the broad-band properties of MKN 231. 



We identify the optically thick absorber with the 
BAL outflows which originate close to the central source. 
The ionized matter, which scatters or reflects the pri- 
mary emission, could be the inner part of this outflow, 
or "shielding gas" (Murray et al. 1995) in agreement 
with the scenario proposed by G02. In both models 
proposed to take into account the AGN emission the 
scattered/reflected components are absorbed. This latter 
absorbing medium could be identified with the starburst 
regions (see e.g. Levenson et al. 2001) or with a different 
line of sight through the BAL wind. 

Beside the flat shape of the observed MKN 231 X-ray 
continuum, the other topic that emerged, since the first 
hard X-ray observations, is the low 2-10 keV luminosity. 
If we assume that the far- infrared emission of MKN 231 
is mainly powered by AGN activity then the intrinsic 
X-ray luminosity should be 10*^ ergs (Risaliti 
et al. 2000). In both models the intrinsic 2-10 keV 
X-ray luminosity can be measured from the absorbed 
component and is ^ l^Jls ^ Although the 

balance between the X-ray and bolometric luminosity is 
derived assuming a standard type 1 AGN, this estimate 
suggests that the contribution of the starburst activity 
to the FIR luminosity could be significantly higher than 
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Table 1. Pn plus PDS data: rest-frame results from the basic reflection-dominated models discussed in Sect. 4.1.1 



AGN Continuum 
MODEL Parameter Value 


Fe line 

Parameter Value 


Nh 1 

A) refl PL+Fe line R 1.3lo i 

iVn / 


E^,„ 6.41+^;^^ keV 
EW < 0.24 keV 

X'/dof 35.61/29 


pa -1 Q9+().;^7 

< / 

B) refl PL+Fe line +scatt PL R 1-Olif*oi 

A7-C 14+1.35 

-''-H ^•-'-^-1.25 


6.42j;E',?g kcV 
EW 0.28i;;'f(; keV 
<T 0.16l°-f° keV 
xVdof 19.6/26 


pa 2.84t°l^ 

ivi o.i9«;ij 

C)ioniz refl+Fe line R 0.94+^^6" 

1 


^^kc 6.40l°:?^ keV 
EW 0.18l!',-l^ kcV 
a O.lOi;^;^^ keV 
X^/dof 102.6/65 


2.951^-1^ 
AT^ 0.24l0;l^ 
D)ioniz refl+Fe line+trans PL R 1.04+o°f5 

AJj 272.6^^;^ 


Bk. 6.39t°:2o keV 
EW 0.18+°o\l keV 
a 0-Wtow keV 
xVdo/ 89.0/61 



For model A and model B we have considered the 2-50 keV data, while for model C and D we have considered the 0.2-50 keV 
data. For model C and D the ionization parameter has been constrained to lie within the range 100-300 org cm/s. 
" intrinsic photon index of the reflected and scattered component; in model D the photon indices of the transmitted and 
reflected component have been tied together; 

absorption of the reflected component in units of 10^^ cm^^; 

absorption of the scattered (transmitted) component in units of lO^^cm^^; 
* unconstrained parameter. 



that of the AGN. 



5.2. The starburst component in MKN 231 

In both AGN models proposed we found a strong soft 

X-ray component which wc describe with thermal emis- 
sion(s) associated to the starburst activity. The derived 
temperature(s) are in agreement with the results obtained 
for the extended emission by G02 using Chandra data 
{kT ~ 0.3 keV and A;r ~ 1.1 keV) and consistent with 
that measured in other well known bona-fide luminous 
starburst galaxies. 

The energetics of this thermal component(s) (L(q 5_2) = 
6 — 9 X 10^^ erg depending on the AGN modelHng) 
is comparable with what has been previously found for 
ULlRGs (Franceschini et al. 2003, Ptak et al. 2003). The 
measured luminosities are in agreement with the one es- 
timated using ROSAT and ASCA data (Iwasawa 1999; 
Turner 1999), but about a factor 4-6 higher than that 
found by G02 using Chandra data (L(o 5-2) = 10^^'^ erg 
s-i). 



For the hard X-ray emission, which we attribute to 
the HMXB, we found that L(2-io) ranges from 6.9 x lO^'^ 
erg s~^ to 1.6 x 10^^ erg s~^ (in the reflection and in the 
scattering-dominated scenario respectively). This spectral 
component is expected in a source like MKN 231 with a 
strong nuclear starburst and has to be taken into account 
in the spectral modelling. We have checked if the X-ray 
luminosity of the cutoff PL component is in agreement 
with the results obtained by Franceschini et al. (2003) for 
a sample of ULIRGs observed with XMM-Newton. These 
observations have shown that for all the starburst dom- 
inated ULIRGs the hard X-ray luminosity of the X-ray 
binaries clearly correlates with the FIR luminosity. 

In Fig. 3 we show the result of the comparisons be- 
tween the X-ray emission due to starburst and to AGN ac- 
tivity and LpiR for the 10 ULIRGs observed with XMM- 
Newton; in these plots we have assumed that MKN 231 
starburst activity contributes for ~ 50% to the Lpju. 

The range of the 2-10 keV HMXB luminosity de- 
rived with the two different AGN models, clearly lo- 
cates MKN 231 in the region occupied by the starburst- 
dominated ULIRGs. This result gives confldence that the 
power-law component originates from the population of 



8 



Braito et al.: XMM-Newton and BeppoSAX observations of MKN 231 



Table 2. Broad-band best fit models and parameters: Model A refers to the reflection-dominated model, while model 
B refers to the scattering-dominated model. Both models include the starburst components. 





Starburst 
Parameter 


Value 




AGN 
Parameter Value 




kT 


u.uu 


+U.U5 
-0.05 


keV 


p 


^•^°-0.11 




/ 










r, 9r|+0.0S 
tt.^U_Q 04 


MODEL A 


Lb/jv 


1.1* 








9fi4- «+l''3.2 
^D4.D_g5 






4.22 


+ 16.9 
-4.22 






6 SS+S JS keV 












EW 


0.20i^;l^ keV 




kTi 


0.35 


+0.09 
-0.05 


keV 


r 


1 OO+0.12 

1.00_Q gg 




kTa 


0.91 


+0.15 
-0.10 


keV 






MODEL B 


Ls/iv 


i.r 








177 5+105.8 










Ska 


6.S9to\l keV 












EW 


0.29l°;J^ keV 



" absorption of the HMXBs component in units of 10 cm^ 

absorption of the reflected/scattered component (model A/model B) in units of lO^^cm"^ 
° absorption of the transmitted component in units of lO^^cm"^; 
the symbol * indicate that the parameter has been kept fixed. 

Table 3. Observed X-ray Fluxes and intrinsic X-ray luminosities derived in the two competing scenarios 

Flux (0.5 10) keV L(0.5 2) keV L(2-10) keV 

10"^^ ergs cm-2 s"^ 10''^ ergs s"^ 10*=^ ergs s"^ 



MODEL TOTAL AGN Starburst AGN Starburst AGN Starburst 



A 7.9 6.3 1.6 33 0.06 20 0.07 

B 8.2 5.0 3.2 2.9 0.09 4.5 0.16 

Model A refers to reflection-dominated model 
Model B refers to scattering-dominated model 



HMXB expected to be present in the circumnuclear star- 
burst region of MKN 231. 

We have also compared the star-formation rate (SFR) pre- 
dicted from the 2-10 X-ray emission with that derived by 
the FIR luminosity. To do that we have taken into ac- 
count that the latter contributes to the 2-10 keV emis- 
sion mainly through the HMXB emission in case of in- 
tense and recent star-forming activity (Franceschini et al. 
2003, Persic et al. in prep.)®. This assumption is consistent 
with a significant contribution of the starburst activity to 
the LpiR. Our result is in fair agreement with the results 
obtained with the modelling of the 1 1000 /xm Spectral 
Energy Distribution by Farrah et al. (2003). 

Considering the heavily obscured PL detected with the 
PDS and the above starburst scenario MKN 231 no longer 
appears as an unusually X-ray faint AGN. In conclusion, 
XMM-Newton and BeppoSAX observations suggest 
that starburst activity strongly contributes to the LpiR 



The relation between HMXB luminosity and SFR hcis been 
recently discussed in the literature (i.e. Grimm et al. 2003; 
Gilfanov et al. 2004 and RanaUi et al. 2003). 



of MKN 231, in agreement with multiwavelength analysis. 



6. Summary 

By combining the XMM-Newton spectral resolution and 

throughput with the high-energy bandpass of BeppoSAX 
we have investigated in deeper detail than previously pos- 
sible the broad-band (0.2-50 keV) X-ray spectral proper- 
ties of MKN 231. 

We confirm that MKN 231 is a complex source whose 
X-ray spectrum has clear signatures of both a highly ob- 
scured AGN and strong starburst activity. 
The detection, in the PDS energy range, of a highly ab- 
sorbed (A^H ~ 2 X 10^^ cm~^) component for the first time 
makes it possible to measure the AGN intrinsic 2-10 keV 
luminosity of MKN 231. This luminosity, which ranges be- 
tween 4.5 X 10^^ ergs s~^ and 2 x 10^^ ergs s~-^, is at least 
an order of magnitude higher than what has been previ- 
ously measured. 

We found that the observed flat 2-10 keV X-ray contin- 
uum of MKN 231 can be explained assuming that below 
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Fig. 2. 0.2-50 keV data coverage. Upper panels: best-fit models for the reflected-dominated (left panel) and scattering- 
dominated scenario (right panel). The spectral components are: a) a highly absorbed PL AGN component; b)left panel: 
a pure reflection AGN component from slightly ionized material or right panel: scattered AGN component; c) a cutoff 
PL component associated with the binaries in the starburst; d) a narrow Gaussian line at 6.39 keV; and e) a thermal 
emission component associated with the starburst. Lower panels: residuals for the different detectors (from up to down: 
pn, PDS, M0S2, MOSl, MEGS). 



10 keV we are seeing only reprocessed radiation through 
reflection and/or scattering. 

The starburst activity significantly contributes to the 
soft (E < 2 keV) X-ray emission (ixh (o.5"2kcV) ~ 3 x 
lO^^ergs s~^), although a relevant contribution due to 
HMXB is also measured in the 2-10 keV energy range 
(-^HMXB (2-iokeV) = 0.7 - 1.6 X lO-^^ergs s-^). 
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